Anisotropic softening of collective charge modes in the vicinity of critical doping in a 

doped Mott insulator 
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Momentum resolved inelastic resonant x-ray scattering is used to map the evolution of charge ex- 
citations over a large range of energies, momenta and doping levels in the electron doped Mott insu- 
lator class Nd2-a;Ce :E Cu04. As the doping induced AFM-SC (antiferromagnetic— superconducting) 
transition is approached, we observe an anisotropic softening of collective charge modes over a large 
energy scale along the V — > (n, 7r)-direction, whereas the modes exhibit broadening (~ 1 eV) with 
relatively little softening along T — > (jr, 0) with respect to the parent Mott insulator (x=0). Our 
study indicates a systematic collapse of the gap consistent with the scenario that the system dopes 
uniformly with electrons even though the softening of the modes involves an unusually large energy 
scale. 
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PACS numbers: 78.70.Ck, 71.20.-b, 74.25.Jb 

The evolution of a strongly correlated material with 
doping— from a Mott insulator to a conducting metal— is 
one of the most intensively studied issues in current con- 
densed matter physics. This fascinating evolution has 
proven to be full of surprises such as the appearance 
of high-T c superconductivity, non-Fermi liquid behavior, 
and nanoscale phase separation [l[ . Mott insulators often 
exhibit phase transitions upon doping, which are signaled 
or hallmarked by the softening of collective charge or spin 
modes. The behavior of spin modes has been investigated 
extensively via neutron scattering Although charge 
excitations near the Brillouin zone (BZ) center can be 
accessed by optical techniques 0, their behavior with 
momentum over the full BZ remains largely unexplored. 
Here, as demonstrated in recent experimental [3, H| an( i 
theoretical studies @, 0] , inelastic x-ray scattering indeed 
provides such a unique opportunity. 

While previous studies have focused largely on either 
the undoped 1-D 8] or 2-D 0, [f| insulators, or the 
hole-doped superconductors 0], Mott insulators can be 
doped with electrons as well. In fact, it appears that 
with electron doping bands in the cuprates evolve in 
a much more straightforward and systematic manner 
[13, EH, H, O 14 1 than with hole doping. However, a di- 



rect study of the doping induced changes in the collective 
excitation spectra of a Mott insulator is lacking. Previ- 
ous x-ray scattering work [l5| focused on the highly elec- 
tron doped superconductor (SC) Ndi.gsCe^o.isCuO^ 
observed intraband and interband excitations and inter- 
preted the excitations in a one-band model with long- 
range hoppings. However, Ref[15] does not report the 
undoped Mott insulator and the way gap modes of the 
Mott insulator are related to the doping where super- 
conductivity sets in. In this Letter, we report a high- 



resolution study of how the collective charge excitations 
of the Mott insulating state (x=0) develop with electron 
doping in approaching the critical AFM-SC transition 
(before reaching the optimal doping) for the first time. 
Our finding, made possible by studying the insulating 
states, is that, as the electron-doping induced AFM-SC 
transition is approached from the x=0 Mott side, the 
system exhibits an anisotropic softening of excitations 
over a large energy scale (~ 5 eV) along the T — > (tt, n)- 
direction, whereas the modes exhibit broadening (~ 1 
eV) with relatively little softening along T — > (tt,0). Our 
results suggest that a multi-orbital Hubbard model is es- 
sential to describe the observed evolution from the x=0 
Mott state in contrast to the 1-band model with long- 
range hoppings proposed in Ref[15] based on the super- 
conductor only. Moreover, our results show that the 
evolution of Mott physics of the electron-doped cuprate 
is dramatically different from that reported in the hole- 
doped cuprates 

The electronic structure of Nd 2 _ a; Ce a: Cu04 (NCCO) 
has been studied by angle-resolved photoemission 
(ARPES) and optical spectroscopies. ARPES studies 
[1 €31 ] find that the electrons dope directly into the bot- 
tom of the upper Hubbard band and yield small Fermi 
(FS) pockets, with a crossover around optimal doping to 
a large FS. Such a scenario, where the magnetic order 
remains commensurate without signs of 'stripe' (charge 
inhomogeneity) or other phase separation, also describes 



magnetization [ 1 3l. 1171 118] and optical data[R 



The experiments were performed at CMC-CAT beam- 
line 9-ID-B at the Advanced Photon Source. Resonant 
inelastic x-ray scattering (RIXS) at copper K edge allows 
a large enough momentum transfer to cover several BZ's. 
The scattered photon energy was measured by using a 
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FIG. 1: Incident energy dependence of electronic excitations: 
(a) Vertical scattering geometry employed with x-ray field 
along [001]. (b) and (c) X-ray energy dependence of inelastic 
excitations for Nd2CuC>4 and La2Cu04. Q- values are chosen 
to minimize large quasielastic background, (d)-(g) Excita- 
tions near the Mott gap are seen to be enhanced near the 
first absorption peak, not only in the undoped insulator, but 
also in the doped system. A few energy-loss curves corre- 
sponding to the data images in (f) and (g) are shown in (h) 
and (i). The red curves show the fits to the quasielastic data. 



diced Ge (733) crystal analyzer, and the intensity was 
recorded by a solid state detector. The overall energy 
resolution was set to about 0.37 eV in order to improve 
count efficiency, enabling us to detect fine momentum- 
dependent details of the spectra. All data are taken at 
room temperature. To avoid possible polarization in- 
duced artifacts when measuring the in-plane anisotropy, 
the sample was mounted with incident polarization di- 
rected along the c-axis (Fig. 1(a)). The data were col- 
lected at several values of the momentum transfer vector 
q = ki — kf in the 2nd BZ along the [100] direction, and 
in the 4th BZ along the [110] direction. The coupling to 
various excitations in the x-ray scattering process, near 
a resonance, depends in general on the energy of the in- 
cident photons. Accordingly, we first examined the de- 
tailed incident energy dependence of the loss spectra at 
several momenta. Representative data sets are presented 
in Fig. 1. It is known that the charge-transfer gap excita- 
tions resonate near the first absorption peak [9j . Similar 
resonance behavior is seen in Nd2CuC>4, which is like that 
in the more extensively studied I^-rrSr^CuC^ 0]. How- 
ever, the scattering intensity at the low-energy branch is 
about an order of magnitude weaker. Our systematic 
investigations summarized in Fig. 1 show that similar 
shape of the resonance profile is seen for the lower en- 
ergy excitations in the electron doped system. Accord- 
ingly, we employed the photon energy of 8.982 keV (this 
energy is above the Is— >3d edge which suppresses the 



crystal-field excitations [20( as it is not the focus of the 
current study). Therefore, in the theoretical simulations 
such incident energy correspondences of the excitation 
spectra are maintained for comparison with the measured 
data. 

Fig. 2 summarizes our RIXS data on NCCO. In the un- 
doped system in panel (a), a broad excitation is observed 
near the zone center around 2 eV with an onset energy 
of rs 1 eV, which is consistent with the charge transfer 
gap found in this compound in optical studies 21 1. It 
is seen to lose intensity as it approaches the zone corner 
along the [ir,ir] as well as the [irfi] directions (Fig. 2(a)). 
Along [7r,7r], it merges at the zone corner with a higher 
excitation band near 5 eV which disperses upward. With 
doping these two high energy branches split at (71", 71") as 
seen in (b) and (c), but only the lower energy branch 
softens significantly (i.e. moves to lower energy), so that 
the overdoped system in (c) displays a large excitation 
gap at the (tt, tt) point from 2 eV to 4 eV. In sharp con- 
trast, the aforementioned 2 eV branch evolves uniformly 
along [tt,0], and softens rapidly near the zone center as it 
tends to close the excitation gap with doping. Notably, 
the zero-loss energy is not accessible due to the presence 
of the strong quasielastic peak (see, e.g., Fig. 1(h)), mak- 
ing it difficult to extract significant data below w 0.4 eV, 
where fitting and subtraction of the quasielastic peak lead 
to uncertainties. Nevertheless, the softening of the low 
energy excitations is evident in the data over the BZ. 

The preceding observations are further highlighted 
through the directly measured data curves presented 
in Fig. 3. Doping dependent changes in the individual 
energy-loss curves are compared in (a). Along [110], dop- 
ing splits the peak in the undoped spectrum around 5 eV 
at (7T, tt) (red uppermost curve on the left side of (a)) into 
two peaks, and the lower of these peaks softens rapidly 
with doping. In contrast, along [100] a monotonic soft- 
ening of the low energy excitations is found. Changes in 
the positions of various spectral features and some of the 
associated leading edges are plotted in (b) as a function 
of momentum and doping. The upper branches around 
5-6 eV are seen to be affected weakly by doping. De- 
spite a significant softening of one of the modes around 
(tt, tt), the closing of the excitation gap is limited to the 
region of the zone center as one approaches the super- 
conducting phase (x = 0.14). These results are in clear 
contrast to the behavior in hole-doped cuprates such as 
the La2- x Sr :E Cu04 or the FBa2Cu307 +( i seriesQ. 

Recent theoretical analyses @, EH of data on x=0.15 
samples suggest two different scenarios. The analysis of 
the data using a 1-band model Ref[15] suggests a no gap 
collapse scenario whereas the analysis of high doping data 
in Ref. [7] suggests a gap collapse scenario. Which of these 
two opposite viewpoints is correct cannot be ascertained 
from the data of Ref. [15] alone. Therefore it is crucial 
to study the evolution of the undoped compound of this 
series with high resolution so that the low-energy gap 
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FIG. 3: (a) Spectral curves (raw data) along [110] and [100] 
for the undoped system (red points) are compared at various 
momenta with the corresponding spectra for the doped sys- 
tem (black and green points) . (b) Positions of peaks (marked 
in Fig.-2(a-c)) based on the center of gravity method 0, [s| of 
various spectral features (filled circles) and some of the asso- 
ciated leading edges (open circles) as a function of momentum 
and doping are plotted. 



FIG. 2: Doping evolution of charge excitations in NCCO. (a) 
x=0 (parent Mott state), (b) x=0.09 (doped AFM), and (c) 
£=0.14 (superconductor). Upper row shows intensity maps in 
the reduced-zone using a color scheme where high intensity is 
denoted by red and low by blue. Black, grey and red dots are 
guides to the eye for the dispersion of low, medium and high 
energy excitation branches, respectively. Lower row gives a 
few representative spectra corresponding to the images in the 
upper row. Absolute momentum values (in units of 7r/a ) are 
shown attached to various spectra and lie in the 4th BZ along 
[110] and the 2nd BZ along [100]. 



excitations can be accessed and model the evolution of 
the spectra from the x=0 Mott side. To interpret the 
present data in terms of RIXS spectra computed within 
the framework of a three-band Hubbard Hamiltonian of 
NCCO, based on Cu d x i_ y i and two O p a orbitals, we ex- 
tend the theoretical framework described in Q by incor- 
porating the doping evolution data made available here 
for the first time. The specific values of the parameters 
used in this work to fit the spectra are : tcuO = 0.85 
eV, too = -0.6 eV, An = -0.3 eV, U p = 5.0 eV; where 
tcuO and too are the Cu-0 and 0-0 nearest neighbor 
(NN) hopping parameters, [n p ] is the average electron 
density on Cu [O], ma the average electron magnetiza- 
tion on Cu, and U [U p ] is the Cu [O] on-site Coulomb 
repulsion. The remaining parameter, crucial for fitting 
the doping evolution of the excitations, the Hubbard U, 
is taken to be 7.45 eV at x=0 with weak doping depen- 
dence: (7=6.69 eV at £=0.09 and £/=6.27eV at 1=0.14, 
so that the effective U decreases by about 16% over this 
doping range, reflecting presumably the effects of screen- 
ing. We note that before proceeding with RIXS compu- 
tations, we determined the chemical potential and the 
magnetization m<j on Cu sites self-consistently at each 
doping level, md values so found are: 0.32 at x=Q; 0.19 



at 2=0.09; and 0.12 at x=0.14. 

Figure 4 (top) shows the calculated RIXS intensity 
maps. The positions of various experimentally observed 
spectral peaks (filled circles) and the leading edge of the 
low energy feature (filled diamonds) are superposed for 
ease of comparison. The high energy peaks involve tran- 
sitions from the nonbonding O and bonding Cu-0 bands 
to unoccupied states in the antibonding band, and fall 
in the same energy range as transitions involving other 
Cu and O orbitals not included in the present 3-band 
model. Therefore, it is appropriate to concentrate on the 
behavior of RIXS peaks within a few electron volts, which 
are associated with the antibonding Cu-0 band, split by 
AFM ordering. In this energy region the theoretically 
predicted changes in the energies of various features as 
a function of momentum and doping are in reasonable 
accord with experimental observations highlighted in the 
discussion of Figs. 2 and 3 above. In particular, the soft- 
ening of low energy peaks is well reproduced and the 
collapse of the gap occurs only near T in the computa- 
tions. The latter effect can be readily understood. The 
RIXS transitions involve both intraband and interband 
contributions. In the absence of a gap, only intraband 
transitions are possible near T which can only exist near 
zero energy transfer close to the Fermi level. Away from 
r, intraband transitions can take place at finite energies, 
but when a gap opens up, interband transitions become 
allowed, even at T (23[. 

Although the computed maps in the top row of Fig. 4 
are shown unbroadened to highlight spectral features, the 
theoretical spectra shown in the bottom panels of Fig. 4 
(blue lines) have been smoothed to mimic the broadening 
of experimental lineshapesj^- The relative intensities 
of the computed peaks are seen to be in agreement with 
our experimental results in lower energies [2J|, excepting 
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FIG. 4: Top row: Doping evolution of RIXS spectra obtained 
within a three-band Hubbard model as a function of energy 
transfer ui at momenta q along F to (tt, tt) and (tt, 0) directions 
for three different dopings x. Positions in uj — q space of vari- 
ous experimentally observed spectral peaks (filled circles) and 
the leading edge of the low energy feature (filled diamonds) 
are superposed. Computed spectra have not been broadened 
to reflect experimental resolution in order to highlight spec- 
tral features. Bottom row. Comparison of theoretical (blue 
lines) and experimental (red lines) spectra. Fhree columns 
refer to the three indicated doping levels x. Each column in- 
cludes four different g-values as shown; experimental spectra 
are at the closest g-value measured. Theoretical spectra have 
been smoothed to reflect experimental broadening 22] . 



the momentum region near the (tt, 7r)-point, where both 
the experimental and theoretical intensities are weak, but 
the computed cross-section is smaller. A similar behav- 
ior is seen in the L^-zSr^CuC^ series [9(. Excitonic 
effects associated with longer range Coulomb coupling 
(intersite-V) beyond that included in our model compu- 
tations could enhance and redistribute spectral weight 
near zone boundaries (25j . The low-energy, gap-edge dis- 
persion relations (E vs. q) ar e in agreement as seen from 
the top row of Fig-4(a-c)[23|]. 

Taken together, these results suggest that the parent 
x=0 Mott insulating state of NCCO is doped uniformly 
by the addition of electrons (within the resolution) and 
that the AFM-gap collapses near the onset of the SC 
phase as seen from the systematics in Fig. 5. In Fig. 5(b) 
we plot the doping dependence of the leading-edge gap 
(edge, green squares) and peak position (peak, white di- 
amonds) of the lowest branch at T. [Since quasielastic 
scattering makes it difficult to obtain experimental data 
right at T, these values are found by extrapolation from 
nearby points, as shown in Fig. 5(c).] Within the exper- 
imental uncertainties the gap collapses in the vicinity of 
the T-point as one approaches the critical doping regime 
in the phase diagram. For comparison, we also show 
in Fig. 5(b) A m = Uirid, the AFM-correlation gap@, 
based on our present calculation (filled circles) and the 




Doping x (%) Doping x (%) q 



FIG. 5: (a) Schematic temperature-doping (T — x) phase 
diagram of NCCO. (b) The leading-edge gap (green squares) 
and peak positions (white diamonds) of the lowest energy ex- 
citation branch in RIXS data at the F-point are compared 
with AFM correlation gap A m from the present calculation 
(filled circles) and from ARPES (triangles) [13] as a function 
of doping, (c) The excitation modes (branches) are fitted and 
extrapolated to the F-point(0,0) based on the branch curva- 
tures, to extract the data points presented in (b). Peak posi- 
tion are shown in the raw data in Fig-2(a-c) and leading-edges 
are from data in Fig-3. 



corresponding ARPES (lol. I2H ] results (triangles), which 
are seen to be in agreement. In general, our analysis 
suggests that the anisotropic softening involving a large 
energy scale observed near the Mott insulating state re- 
quires going beyond the one band model. 

In summary, we have utilized the unique momentum 
resolution of x-ray scattering to map the evolution of 
particle-hole excitations in the electron doped cuprate 
from its parent Mott state (x=0) to the doping on-set 
for superconductivity. We observe a nearly degenerate 
charge excitation mode near the (tt, tt) point around 5 
eV in the Mott insulator, which splits on doping away 
from half-filling, with its lower energy branch softening 
anisotropically in approaching the AFM-SC critical dop- 
ing. In contrast, the response near the (7r,0) wavevector 
exhibits damping with relatively little softening. Our re- 
sults indicate a systematic collapse of the gap consistent 
with the scenario of uniform doping when electrons are 
added into the parent Mott state which is in clear con- 
trast to what is reported in hole-doped Mott insulators 
where charge inhomogeneity is commonly observed. The 
unusual softening of the modes from the Mott state likely 
contributes to many unusual behaviors of cuprates. 
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